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We introduce the notion of general K, = — 1 Friedman-Lemaitre (compact) cosmologies and 

the notion of averaged evolution by means of an averaging map. We then analyze the 

Friedman-Lemaitre equations and the role of gravitational energy on the universe evolution. 

We distinguish two asymptotic behaviors: radiative and mass gap. We discuss the averaging 

problem in cosmology for them through precise definitions. We then describe in quantitative 

detail the radiative case, stressing on precise estimations on the evolution of the gravitational 

energy and its effect in the universe's deceleration. Also in the radiative case we present a 

smoothing property which tells that the long time H x H stability of the flat JC = — 1 

FL models implies W +1 x H l stability independently of how big the initial state was in 

I H l+1 x H 1 , i.e. there is long time smoothing of the space timer] Finally we discuss the 

5— ( 

(3JT) existence of initial "big-bang" states of large gravitational energy, showing that there is no 

mathematical restriction to assume it to be low at the beginning of time. 

(N 

"The word smoothing here is referred to the decay toward zero of the space time Bel-Robinson 
C ~ ) curvatures (and therefore of the derivatives), and not to a gain in Sobolev regularity as in usual PDE 

terminology. 

o 

ON 

o 

1 Introduction 

O 

m ^ An implicit assumption of the Friedman-Lemaitre cosmologies as models of the actual universe is 

that, because matter distribution at large scales (visible or not) appears to be "to a good extent" 

H 

homogeneous and isotropic, the large scale evolution of the universe should be modeled as driven 
"to a good extent" by an exactly homogeneous and isotropic material distribution. The assumption, 
now known as the averaging problem in cosmology, needs quantitative approval or disproval (see [1]). 
Phrasing the problem in a question one asks: is the large scale evolution affected by the small scale 
structure?. The reason of the difficulty lies evidently in the nonlinearity of the Einstein equation. An 
averaged source of matter doesn't give rise necessarily to the average of the original solution. We will 
discuss this and other issues from the perspective of general cosmological models, i.e. the study of 
arbitrary solutions of the Einstein equation in the Hubble gauge (constant mean curvature (CMC) 
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gauge) provided with a set of Friedman-Lemaitre equations giving the cosmological interpretation to 
the framework. 

The standard JC — — 1 FL cosmology describes the universe history by the evolution of the en- 
ergy and pressure densities of the different type of matter present. Starting from a "big-bang" where 
the densities and the space-time curvature blow up, the universe evolution is described as eternally 
expanding, with decaying densities and space-time curvature at a particular pace according to their 
matter type. Such a description is analytically possible due to the homogeneity and isotropy of the 
space which reduce the Einstein equations into a set of ordinary differential equations, the so called 
Friedman-Lemaitre equations. We will deal here with compact cosmologies, i.e. space-times with 
compact Cauchy surfaces of hyperbolic type. When speaking about homogeneity and isotropy of a 
compact cosmology we will refer to those properties in the universal cover solution. In its formal terms 

2 

the geometric structure of the space-time is described by a metric of the form g = —dr 2 + a 2 (r) /V^gn 
on a four-manifold R x £ where £ is a compact hyperbolic manifold i.e. a manifold admitting a metric 
of constant negative sectional curvature and where Vr is the volume of T, endowed with the unique 
hyperbolic metric (the one with sectional curvature equal to negative one). If the densities of energy 
and pressure of the material fields are p(r) and p(r) the FL equations are 

2 8nGp JCVI 

n = —3 r (1) 

a" _ 4ttG( P +3p) 

a ~ 3 (2j 

where H = a'/a is the Hubble parameter and G is the gravitational constant. These equations must 
be complemented with an equation of state p(p) . An obvious observation about these models is that 
they do not have any pure gravitational degree of freedom besides the gravitational field generated 
by the matter present. This fact is seen by making p = p = and observing that in that case the 
solutions are flat. We call these flat solutions flat cones as they can be obtained as quotients of a future 
light cone in Minkowski space-time. For non homogeneous and isotropic solutions there is no way to 
define which part of the gravitational field is generated and which part is free, as those properties (if 
anything) would be potentially defined only in special solutions or in asymptotic regimes. There are 
simply two fields interacting, gravitation and matter. In this sense, the gravitational field adds a new 
degree of freedom to general cosmological models which needs to be quantitatively described. 

We have found that a satisfactory way to analyze arbitrary solutions to the Einstein equations 
on the light of the questions raised by cosmology and those raised by the FL models themselves, is 
to introduce the notion of general cosmological model: an arbitrary solution to the Einstein equations 
in the Hubble gauge, provided with a set of Friedman-Lemaitre equations giving its interpretative 
cosmological meaning. Unlike the FL models where the FL equations are enough to describe the 
evolution, in general cosmological models one must rely on the full Einstein equations to predict 
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the behavior of the terms involved in the general FL equations and therefore interpret the solutions 
in cosmological terms. One purpose of the article is to start a rigorous analysis of the general FL 
equations using the full Einstein equations. 

An arbitrary solution g on a space-time manifold Rx£ where £ is a compact hyperbolic manifold, 
is in the Hubble gauge if the mean extrinsic curvature k of the equal time Cauchy surfaces is constant. 
The foliation R x £ is called the CMC foliation. It is well know to be unique, intrinsically defined, an 
with the mean curvature k varying monotonically on it, in particular k or Ti (as we will se the Hubble 
parameter is TL = ^p) can be taken as a time variable. It is important to remark that unlike other 
gauges, the Hubble gauge is intrinsic, i.e. it is implicitly given by the solution. Let's write the metric 
as 

g = -(N) 2 dk 2 + X* x dk + dk^X* +g, (3) 

where N is the lapse function, X the shift vector and g is a spatial three-metric on E. To write the 
general FL equations one defines the radius a(k) at the time k as a(k) — V(k)s and the proper time 
r(k) at the time k through (see [5] for a related approach) 

dr = jNdv g 

dk v ■ { ' 



With these definitions the FL equations (deduced from the Einstein equations, see subsection 3.1 1 are 



! pv , pvr ,. ^2 SyMlttVg ^N{l^Gp+\K\ 2 )dv g 

1. iirst iL equation : n = — — , (5) 

„ Am ,. a" -^Kf{^G{p + Zp) + \k\ 2 )dv g 

2. Second FL equation : — = — . (6) 

a 3V 

Where Af = (bar denotes volume-average) and has average equal to one. The derivatives, denoted 
with a prime are proper time derivatives, i.e. ' = -4-. K is the traceless part of the second fundamental 
form K. Compared with the second FL equation [2] in a perfect FL cosmology we observe the appear- 
ance of the weight term Af which inexorably couples matter to gravitation and a purely gravitational 
term of \K\ 2 which is essential and represents the additional gravitational degree of freedom mentioned 
before. A particular solution is a FL model iff K = and JV = 1. 

On the light of general cosmological models, a fundamental question is to quantify the evolution 
of the different terms that appear in the FL equations. It is important to realize that the ultimate 
goal would be to understand the FL equations for solutions which are realistic at small scales, i.e. at 
the natural scale of the flow. This is a difficult problem, however we will argue that we can have an 
starting point if precise assumptions are made. Namely, in subsection |3.5| we will introduce assumption 
(C), a precise quantitative hypothesis on the behavior of arbitrary solutions at late times, from which 
we will make explicit estimations of the different terms involved in the FL equations. Assumption 
(C) is a close relative of the weak cosmic censorship conjecture of Penrose, conjecture stated in an 
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asymptotically flat context. In rough terms, assumption (C) precisely describes a family of solutions 
and divides it in two classes: radiation and mass gap. A radiative solution is an ideal solution in 
which no sort of compact object emerges along evolution, i.e. universes filled only with radiation. We 
will study this case in detail, although only for gravitational radiation. The technique may be applied 
to other radiative contexts as well. In this case the gravitational field can safely be isolated from 
the rest, and one can safely interpret j 67rG — Pg — Pg a s the effective energy and pressure densities 
of gravitational radiation. These densities are quantitatively studied along with the decay of Af to 
one. The estimates are given in Theorem 1 (see statement below) which in addition give estimates 
on the Bel- Robinson energies Qi. Altogether Theorem 1 provides a detailed structure of the radiative 
solutions and is one of the main result of this article. 

Theorem 1 (Expansive smoothing and energy estimates). Let £ be a compact and rigid hyperbolic 
manifold. There is an e > such that the Einstein CMC flow of a cosmologically scaled initial state 
(i.e. with Ti = 1) (g,K) with V — V, n / < e and £\ < e has the following long time properties (take 
t = 

1. The limit limt^oot^Qa is finite and greater than zero. 

2. there are rij > such that lim^oo jj^jyr-Qi ^ 00 f or * — !• 

3. for given 7 > 0, J t °° ^^^ du > CH 2 +^). 

4- \K\ 2 < Ct~ A pointwise (not volume averaged). 

In particular the cosmologically scaled flow of a H l x H % state (for any i > X) as in the hypothesis 
above converges in H l x H 1 ^ 1 to the canonical flat cone state (g,K) — {gH,—gu)- 

Theorem 1 is in PDE terminology a small data statement. The small data condition is stated as 
saying that the reduced volume V = H 3 V is e-close to its infimum and the first Bel- Robinson energy 
E\ e-close to zero. These two conditions can be seen to be equivalent [3] to the statement that the 
initial data (g,K) is close in the Sobolev space H 3 x H 2 to the flat cone state (gn, ~9h) where gn 
is the unique hyperbolic metric (up to diffeomorphism) . A hyperbolic manifold is called rigid if it 
doesn't admit traceless Codazzy tensors (see [B] for a discussion). The topological condition of rigidity 
is important to get the precise estimates above. It is possible to get estimates in the non rigid case 
but they are different, in particular those on the gravitational energy. The importance of rigidity is 
that it allows the control of the H 2 norm of harmonic metrics with respect to the hyperbolic metric 
(spatial gauge) only by their Ricci tensor. 

The estimates in theorem [l] are compatible with what one would expect is a radiative behavior. 
According to the standard FL models an exact radiative behavior would imply a pointwise decay on 



4 



the gravitational energy density of the form 

(7) 

l6irG t 4 1 ' 

The estimate in items 3 and 4 in theorem [TJ says that in some averaged sense the global gravitational 
energy decays with a rate between the radiative t~ x and the faster t~ 2 . It would be interesting to 
improve (if possible) the estimate from below in item 3. 

In rough terms the mass gap solutions can be described as those for which after a sufficiently 
long time there appear a finite set of isolated stationary solutions separating from each other and 
with radiation in between. This qualitative description is made quantitative in assumption (C). We 
analyze the averaging problem for these mass gap-solutions. A convenient setup for the analysis is to 
define the notion of averaged space, a Lorentzian manifold constructed out the averaged parameters 
a(k) and r(fc) of the original solution. The averaging problem can be stated as asking to which extent 
the averaged space is close to a FL model. A remarkable consequence of applying assumption (C) is 
that the second FL equation is estimated as 

a" = 4itG(Madm + P + Pg + HP + Pg)) h ^- (3+6 )^ ^ 

where Madm is the volume average of the ADM masses of the emerging stationary solutions, and 
p, pQ , p, pa are the volume averages of the densities of energy and pressure of material and gravitational 
radiation respectively, filling the space in between. However one must remark that despite all the 
satisfactory equation [8] may look, it is based on an idealized assumption and on its apriori estimates 
which so far need to be justified. Also the quantitative description they provide is only asymptotically 
in time, and not throughout the full evolution. 

A quantity underlying all the averaging formalism is the so called reduced volume [5], defined 
above as V = H 3 V. It decreases monotonically, and is bounded below by the topological invariant 
Vh ■ It has been used in [3] to show the long time geometrization of the Einstein flow under curvature 
bounds. Here it is manifested throughout the article in different forms. Its monotonicity is shown to 
be equivalent to the universes deceleration and is used to get the estimate in item 3 in Theorem 1. 
We will introduce an use an equivalent quantity that we will call the global CMC energy defined as 

EcMC = ^ R {V-V in} ). (9) 

Rather remarkably, the CMC energy is shown to express the full ADM energy of the time-asymptotic 
evolution only in terms of the total volume, the Hubble parameter and the topological invariant Vi n f. 

The contents and sections are organized as follows. In section [2] we introduce the main equations 
for the Einstein-CMC flow as well as Bel-Robinson energies and their main formulae. In section [3] we 
introduce the averaged cosmological parameters and the Friedman-Lemaitre equations. The treatment 
has no restriction on the sort of matter. We introduce the Newtonian gravitational potential 0, its 
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Poisson-like equation and reformulate the FL equation with it in subsection |3.2| As it turns out 
the Newtonian potential is the main field to estimate when the purpose is to estimate the universe 



deceleration and the Hubble parameter as a function of red shift z. In subsection 3.3 we introduce 
the CMC global energy and relate it in subsection |3.4| with the ADM energy in the weak field limit, 
analysis extended in subsection 3.5 to arbitrary solutions under assumption (C). In subsection 3.6 we 



discuss the averaging problem on the light of assumption (C) for the mass gap regime. We will use the 
CMC energy to estimate the gravitational energy in section [4] Also in section [4] we prove the main 
estimates of theorem [l] The technique may be thought as estimating the gravitational field through 
a Taylor expansion (in time) of the zero-order Bel-Robinson tensor and is a natural extension of the 
analysis in ([3]). In section [5] we construct "big-bang" states of high gravitational energy showing that 
there is no mathematical reason to assume a low gravitational energy at the initial "big-bang" state. 
The dynamics of those states even in short times is a completely open problem, in particular it is not 
know whether the initial rate of expansion with respect to proper time is of matter, radiation or of a 
type like non of both. In section [6] we give an account of the main points of the article. 



2 The CMC flow equations and the Bel-Robinson energies. 

2.1 The CMC flow. 

In this section we consider the formal setup of the Einstein CMC flow equations. A detailed account 
can be found in [3]. Consider £ a compact hyperbolic three-manifold. A cosmological solution to the 
Einstein equations with compact Cauchy surface S is formally a Lorentz metric g on a four manifold 
of the form / x S (where / is a interval) and where the equal time hypersufaces E t are space- like i.e. 
the induced metric is Riemannian. If the mean extrinsic curvature (fc = tr g K) is constant on each 
slice of the foliation {S t } then we say that the cosmological solution is in the (temporal) CMC-gauge. 
When the spatial topology is a hyperbolic manifold the mean curvature k cannot be zero (due to the 
energy constraint and the fact that E doesn't accept metrics of non-negative scalar curvature) and it 
can be proved to be strictly monotonic over a unique and connected interval. For a three manifold 
of hyperbolic type in particular it is conjectured that the CMC foliation has a range of k equal to 
/ = (—oo,0), i.e. from a "big-bang" when Ti. — * oo towards an infinitely expanding universe when 
H — ► 0. Say dt — NT + X where T is the unit normal to the slices and t = k. Write the four metric 
as 

g= -N 2 dt 2 + X* ®dt + dt® X* +g, (10) 
where g is the spatial three dimensional metric. N is called the lapse and measures the rate of proper 
time with k (locally). X is called the shift vector field and can be chosen freely but compatible with 
the regularity. For a discussion of the initial value formulation in the CMC gauge we refer the reader 
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to [3]. We call the path (g, N, X)(k) the CMC flow. A CMC state is a pair position-normal velocity 
(g, K) (where K is the second fundamental form and is equal to K = — ^Cxg) with k = tr g K constant. 
Thus the CMC flow gives rise to a flow of position and velocities (g,K)(k). With this notation the 
Einstein equations 

Rice - ^Rg = 8ttGT, (11) 
can be seen as the CMC flow equations (taking t = k) 

1. Hamilton- Jacobi equations 

g 1 = -2NK + C x g, (12) 
K 1 = -S7 2 N + N(Ricc + kK-2KoK)+ £ X K - 8nGN(T - ^J^g), (13) 

2. Constraint equations (energy and momentum respectively) 

R - \K\ 2 + k 2 = 16ttGp, (14) 
V.K = SttGJ, (15) 

3. Lapse equation (deduced from equations above) 

- AN + (4ttG(p+3 P ) + \K\ 2 )N = I. (16) 



The T-term in the right hand side of equation 13 must be thought to be restricted to E. Also as usual 
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p = T(T, T), J — T(T, .) and p = ( Tab ^- 9 — I is the average of the principal pressures. In equation 
V.if = \7 a K a b is the divergence and in equation 13 it is (K o K ) a b = K ac K c b . Finally the speed of 
light was taken to be c = 1 . 



2.2 The Bel- Robinson energy and the space time curvature. 

We will measure the L 2 norm of the space time curvature relative to the CMC gauge. We will also 
need to measure the L 2 norm of their time derivatives relative to the normal direction to the CMC 
foliation. There is a remarkable way to introduce them and it is by means of Weyl fields. Although 
we won't discuss Weyl fields in detail as there are very accurate references on the subject ([7], [6]), we 
will mention the most used properties here and briefly elaborate on their conceptual importance as 
variables controlling the gravitational field. 

Definition 1 A Weyl field is a traceless (4, 0) space time tensor satisfying the symmetries of the 
curvature tensor Rm. We will denote them by ^V a bcd 0T simply W. 
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The Ricmann tensor of a vacuum solution to the Einstein equations is a Weyl field that we will 
denote as Rm = Wo. Let T be the normal field (future pointing) to the CMC foliation. Then 
VyWo = Wj are Weyl fields. Together with the volume radius ([5]) and the L 2 norm of the second 
fundamental form K they are an important set of variables that control the gravitational field (i.e. 
the metric g relative to the foliation) see [3 . A central advantage for taking them as variables is that 
they enjoy remarkable algebraic properties that simplifies the space time algebra considerably. We 
discuss the main formulae below. Given a Weyl tensor W define the left and right duals *~W a bcd — 
^abi m W lm cd and W* abcd = W ab lm ±e lmcd . Both are Weyl tensors, *W = W* and *(*W) = -W. 
Define the current J(W) and its dual J*(W) as 

V a W abcd = J abc (W), (17) 

V a W* ahcd = J Q * bc (W). (18) 

For the Riemann tensor in a vacuum solution to the Einstein equation we have J — J* = due to the 
Bianchi equations. This is a central fact that will be of fundamental importance latter. We also have 

V [a W bc]rfe = ±e fabc J*/ e (W), (19) 

V [a W* c] de=ie /abc 4 e (W). (20) 

The L 2 norm with respect to the foliation will be defined through the Bel-Robinson tensor. Given a 
Weyl field W its Bel- Robinson tensor is 

Q abcd (w) = w alcm w b i r + w: icro wv,r- (2i) 

It is symmetric and tracclcss in all pair of indices and for any pair of timelike vectors T\ and T-2, 
Q(Tx, Ti, T2, T2) is positive if W ^ ([7 ). In particular we define the L 2 norm of W with respect to 
the foliation as Q(T,T,T,T). It is seen to be the L 2 norm of the electric and magnetic fields of W 
defined through 

E ab (W) = W acbd T c T d , (22) 

B ab (W) =* W acbd T c T d . (23) 

i.e. Q(T,T,T,T) = \E\ 2 + \B\ 2 . They are symmetric, traceless and null on the T direction. For the 
Riemann tensor in particular we have 

Eab(W a ) = Ricc ab + kK ab - Ka c K c b , (24) 

and 

ejB lc (W ) = V a K bc - V b K ac . (25) 
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The following formulae provide the components of a Weyl field with respect to the CMC foliation in 
terms of the electric and magnetic fields k, I are spatial indices) 

W iifeT = -e t / n B mk (W), *W ijkT = e ZJ m E mk (W), (26) 

W tjkl = e ljm e kln E mn (W), *W ij]d = e lTm e kln B mn (W). (27) 

The divergence formula 

V a Q(W) abcd = W b m d " J(W) mc „ + W ™ " J(W) mdn + (28) 

+* wr;j*(W)™c„ +* W? c n J*(W) mcn , (29) 

and therefore 

V«Q(W) qTTT = 2EV(W) J{W) iTj + 2B i U*(W) iTj (30) 
gives the Gauss equation 

d kQi T ^T,T)dv g = _ ^ 2NEij(yV j j (w) . T . + 2NBVJ*(W) iTj + 3NQ abTT U ab dv g , (31) 

where H ab =V a Tb is the deformation tensor and plays a fundamental role in the tensor algebra. In 
terms of the electric and magnetic fields the components of Q a bTT are written as 

Q iTTT = 2(E A B) t , (32) 

Q ljTT = -{E x E)ij — (By. S) y + |(|^| 2 + \B\ 2 ) 9l] . (33) 

Controlling J and J* in L 2 and II in ff 2 is enough to control the L 2 norm of the Weyl field. The 
following formulas are essential when it comes to get Sobolev estimates of the Weyl field 

divE(W) a = (KAB(W)) a + J TaT (W), (34) 

divB(W) a = -(KAE(W)) + J£ aT (W), (35) 

curlB ab (W) = E(V T W) ab + §(£(W) x K) ab - \kE ab (W) + J aTb (W), (36) 

curlE ab (W) = B(V T W) ab + §(B(W) x K) ab - \kB(W) ab + J* aTb (W), (37) 

where the operations A, x are defined as 

(A x B) ab = e a cd e b ef A ce B df + \(A- B)g ab - l -(trA)(trB)g ab , (38) 

(A A B) a = e a bc A b d B dc . (39) 

The equations 34p7 above are an example of the so called elliptic Hodge systems ([7]). In particular 
under basic regularity of the background metric they make possible to get elliptic estimates. 
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2.3 Scaling. 



Scaling is the operation allowing us to speak like "looking the system at a particular scale" . It is 
a different operation of than coordinate scaling, as scaling a solution does changes the solution but 
scaling coordinate systems doesn't. Both transformations are however important when used simulta- 
neously. 

Definition 2 Given a solution g to the Einstein equations, we call A 2 g the solution g at the scale of 
j and we call A the scale factor. 

We say that a CMC state (g, K) is cosmologically scaled (or normalized) if k — — 3 or the same H = 1 
as we will see the Hubble parameter TL is equal to Given a state (g,K) that gives rise to a 

global solution g we can scale it as ^-g to transform the original state (g, K)(k) into a cosmologically 
normalized state {^-g,=^K). Therefore a state (g,K) has a cosmological scale of ^ = ^. Say 
(g,K)(k) is a CMC state, and say U is some space time tensor constructed out of g that we are 
looking at the fc-slice. The corresponding values of U on the same slice when we cosmologically scale 
the state {g, K) will be denoted with a tilde (cither above or next to it) say U or ?7~. Thus g = ^-g 
and K = ^K. In a CMC flow (g, K){k) we can cosmologically scale the solution g at every k getting 
thus a flow of normalized states (g, K)(k). In the flat cone case the cosmologically scaled flow is just 
(gn, —9ii)(k) and what we will call stability of the flat cone will be the stability of the cosmologically 
scaled solutions. In general a space time tensor will scale as X S U for some weight s, therefore U will 
be just U = {=£) S U. We will indistinctly use -jp or H as the scale factor A. The following table shows 
how some main tensors transform when g — > A 2 g. 



g 


A 2 g 


9 


A 2 . 9 


K 


XK 


k 


k 
A 


N 


A 2 iV 


<f> 
Wi 


A-*+ 2 W 4 


Qi 


A-( 2l+1 )Q, 



4> is the Newtonian potential defined below. 
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3 Averaged evolution. 



3.1 Averaged cosmological parameters and the averaging map. 

We define the geometric parameters, a(k) (universe's radius), r(fc) (proper time), and H(k) (Hub- 
ble parameter) in volume average. All those parameters reduce to the standard FL parameters when 
the solution is homogeneous and isotropic. 

Definition 3 Given an arbitrary CMC solution we define the universe 's radius at an instant of time 
i 

k as a{k) = V3 k *. . The volume-averaged proper time r(fc) is defined through 

Tk = ^v-*- (40) 

Recalling that in the FL models the Hubble parameter is defined asH=-^ we compute 

1 d V^ 1 dVi dk 1 1 2, f Arl , , V -k 

H = — r— i— = -^—r^-r = ~^t,V^^{ / -Nkdv„)- r — = — . (41 

Vh dr vl dk dr ^ 3 V E 9 'J s Ndv g 3 v ; 

Thus in arbitrary solutions TH. = = M. This expression is valid also locally in the following sense: 
define the cube of the local radius as the volume clement dv g (k), then the local Hubble parameter 
is one third the logarithmic derivative of the volume clement with respect to the proper time in the 
normal direction to the CMC slice k. A direct computation gives for the local Hubble parameter 

n_/ _ 1 dvg _ ^fe 
3dv g dr 3 

The Friedman-Lemaitre equations take the form 

1. First FL equation : H 2 = ~ 6y + — &y , (42) 

2. Second FL equation : — = — . (43) 

CI ov 

Where J\f = ~ (bar denotes volume-average) and has average equal to one. The derivatives, denoted 
with a prime are proper time derivatives, i.e. ' = The first FL equation is just the volume average 
of the energy constraint 

\&tt P = R-\K\ 2 + \k 2 . (44) 

o 

Observe that to make it look closer to the second FL equation, we can multiply the energy constraint 
before integrating by Af and integrate thereafter to get 

J s MRdv g J s M(l67rGp+\K\ 2 )dv g 

To obtain the second FL equation we observe that 

(")'=--(-) 2 = --^ (46) 
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and 



( a' y= dH = ldk = V 

{ a ' dr 3dr 3 / E Ndv g ' { 1 



On the other hand integrating the Lapse equation 16 we get 



/ N(iirG(p + 3p) + \K\ 2 )dv g = V - 3H 2 [ Ndv g . (48) 

Equations [46] [47] [48] together give equation [43] 

Let's restate the standard JC = — 1 FL models on the light of the description given above for 
arbitrary solutions. If the solution is g = —dr 2 + a(r) 2 gH on a manifold IxE then a(r) = (y-)^ 
where Vh is the volume of £ with the hyperbolic metric gn and V is the volume with the metric 
a(r) 2 gH- Our choice of radius for arbitrary solutions has been instead a(r) = 1 we will make this 
choice in equations [49] and [50] below . We recall too that in the standard FL models the energy density 
and pressures are a function only of r and for that reason they coincide with their volume averages. 
Taking these facts into account the standard FL equations are 

1. H 2 = S^*"G P )dv g kvI 
6V a 2 



a" _ -/ s (47rG(p+3j>)d« 9 



a W 



Observe that in the FL equation 



(50) 
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instead of the curvature term — JCVjj jar we have the term 



Jg Rdv 9 = / Js-Fftfog x 1 
6V 1 '" 2 



(^TTr^)^^ (51) 



where the first factor in the last term of the previous equations is scale invariant and therefore equal 

2 

to Vjj for any metric scaled from the hyperbolic metric (so is close to it for any metric scaled from a 
metric close to a hyperbolic metric). 

In order to establish a mathematical definition of the averaging problem in cosmology we define 
the averaging map from arbitrary CMC solutions into Lorentzian manifolds in the following way. 

Definition 4 Given an arbitrary CMC solution g on M x S with S a compact hyperbolic manifold 
define the volume-averaged solution as the Lorentzian space (K x E, .4(g) ) with -4(g) = g = —dr + 
a ^ T l 9h, where r and air) are the averaged proper time and radius as given in definition\A gn is the 
unique (up to diffeomorphism) hyperbolic metric that £ accepts. 

It is essential in the definition above that, due to Mostow's rigidity, there is one hyperbolic metric up 
to diffeomorphism in a given hyperbolic manifold. That makes the definition of .4(g) unambiguous. 

In rough terms the averaging problem for arbitrary solutions can be stated as to whether the 
averaged space .4(g) is "asymptotically in time close" to an exact K. = — 1 FL solution with the 
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"averaged energy density and pressures" "asymptotically in time close" to the energy density and 
pressures of the exact FL model. One may also replace "asymptotically in time close" simply by 
"close" all along evolution. Physically that would be a more adequate question to ask. This definition 
however faces various indcfinitcncss, we comment on them below. 

i) The first is to give a precise meaning to "averaged energy and pressures" for arbitrary solutions. 
We can safely say what they are for the material fields, as material fields posses densities of pressures 
and energy, but it is not known what they are for the gravitational field, and presumably they can 
not be isolated as densities. The old question on how to define the gravitational energy which shows 
up all through General Relativity is present also here. A consensual definition of energy is the total 
ADM energy, a global term comprising the energetic content of a global system. Despite all the 
satisfactory the expression is, it is defined in asymptotically flat space-times and not in the context 
of cosmological solutions. We will argue in subsection |3.6| on the validity of the averaging problem , 
at least asymptotically in time, if it is assumed a compact and extended relative of the weak cosmic 
censorship conjecture of Penrose, conjecture stated for asymptotically flat space-times. Indeed we will 
analyze the averaging problem under the assumption that, under a particular model for matter at 
natural scales (the small structure), it happens that, generically, cosmological solutions evolve into a 
finite set (however large it may be) of asymptotically flat stationary solutions separating from each 
other, with gravitational radiation in between and if in addition we compute the "averaged energy 
density" as the volume-average of the ADM energies of the stationary solutions plus the volume- 
average energy of the gravitational radiation in between. Both terms, as we shall see, can safely be 
computed. We will call the assumption above assumption (C). The extent as to whether this idealized 
assumption would be applicable to the actual universe in which we live at present times is not under 
consideration here. However I would like to point out one aspect that immediately jumps out and that 
it would have to be addressed with care. Assuming that galaxies conform the individual stationary 
solutions, there is the issue to establish, due to the large dark halos extended over diameters many 
times their visible diameters, where (if somewhere) and how far the individual galaxies (including their 
halos) become asymptotically flat. This lack of asymptotic flatness on large neighborhoods around the 
visible galaxy is manifested in the well known flat rotation curves of stars with large orbital radius. 

ii) A second problem in the rough definition of the averaging problem given above is to specify 
the equation of state of the exact FL solution from the original solution at natural scales. On the 
light of assumption (C) there are two situations possible, a radiative regime, of universes filled only 
with gravitational radiation, and a massive regime, of universes where in addition to radiation there 
are massive compact object (the stationary solutions). Both regimes, that we will call radiation and 
mass gap respectively, deserve different technical analysis. We will discuss the radiation regime in 
rigor and detail below in section [3] The analysis of the mass gap regime is done in subsection |3.5[ 
Although rigorously deduced from the assumption (C), it lacks a precise determination on the decay 
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of the radiation term. We will return to this point later. 

iii) A third problem is to define in a quantitative manner the notion of "closeness" between the 
averaged space and the exact FL solution. Precisely, we have to specify the scale in which the solutions 
are compared and a law for the asymptotic relation between them. 

3.2 The Friedman-Lemaitre equations and the Newtonian potential. 

A remarkable fact about the averaging formalism is that the second FL equation can be written 
only in terms of the volume-average of the Newtonian potential <p and consequently a(r), 7i(r) and 
z(t) are determined only from <f>. 

Definition 5 Define the Newtonian potential tf> as (f> = 1. It satisfies the Poisson equation 

(Lapse equation) 

A0 = (4ttG(p + 3p) + \K\ 2 ) + {AirG{p + 3p) + \K\ 2 + 3i? 2 )0, (52) 

or making e = 4irG(p + 3p) + \K\ 2 

A(j) = e+ (e + 3H 2 )(/>. (53) 

From the Maximum principle it is seen that — 1 < <f> < 0. Observe too that <fi is an absolute potential, 
i.e. there is no ambiguity in the level of energy in its definition (as can be deduced from the unicity 



of solutions in equation 52) and observe also that it is scale invariant. As defined here the Newtonian 



potential of course coincides with the usual Newtonian potential in the weak field Newtonian regime 



(when p w and K rs 0). Compare also equation 52 with the usual Poisson equation in Newtonian 
dynamics 

A(f> = AirGp. (54) 

Equation [52] is fundamental to understand the dynamics of the gravitational field in general and its 
analysis extracts among other things the time at which Newtonian dynamics appears, i.e. when is that 
gravitation gets ruled by classical Newtonian potentials at large scales. A straightforward calculation 
gives 

-=H 2 -^, (55) 



or 



W _ -1 
H 2 ~ IT 



, , ( 56 ) 
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where cj> is the volume-average of <f>. This equation can be used to get an equation for TL as a function 

i 

of red shift 1 + z — V3 1 (V is the present volume and Viz) is the volume at the corresponding red 

V(z) 8 

shift). The relation is 

d ^ H 1 (57) 



d\n(l + z) 1 + 0' 
One also obtains 
d\n{l + z) 

dr = (58) 
Of course an estimation of </> as a function of r, z or TL is needed to make use of the equations above. 

3.3 The CMC energy. 

We would like to define a formal quantity on CMC states on a compact manifold £ analogous to 



the total ADM mass of asymptotically flat space-times. Restate the first FL equation 49 in the form 



l-( V ^) l =« m , (59) 

where we have defined Vinf as the absolute infimum of the reduced volume V = TL 3 V(g, K) among the 
set of all CMC states (g,K). It is known [5], [5] that if S is hyperbolic Vi„/ = Vjj- Q m is defined as 
usual as TO = 8 3 ^ 2 P . Thus the density of mass p and the Hubble parameter Ji determine the deviation 
of the reduced volume from its absolute infimum. If ~0we get in particular the approximation 

WW ^- V ^ (60) 
This remarkable equation, expresses the total mass M. in terms only of Tt, G, the total volume V and 
the topological invariant Vh- As we shall sec in section [3.5| it holds too, asymptotically in time, for 



general models under assumption (C) . Inspired on it and equation 60 we define the total CMC energy 
as 

Definition 6 Define the CMC global energy as 

E C MC = ^ R (V-V inf ). (61) 

3.4 The ADM limit of the CMC energy: radiation. 

Recall that the Hessian of the ADM energy around the flat Minkowski space-time state g = gE 
and K = (gE is the euclidean metric) is (see for instance [8]) 

8wG5 ( - 2 ' ) E A dm=1 I \\7g' TT \ 2 dv+ [ \K' TT \ 2 dv + %TtG ( 5 {2) pdv, (62) 
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where TT means transverse-traceless with respect to the flat metric The Hessian of the reduced 



volume V was calculated in [9 . We include below a calculation of the Hessian of the CMC energy 61 



based on their analysis for the sake of completeness and clarity. The Hessian of the CMC energy in 
the limit when k — ► is locally the same as equation [62] the precise expression is 

SttGS^Ecmc = [ \K' TT \ 2 dv g + l - f \Vg' TT \ 2 dv g - ~ f \g' TT \ 2 dv g + 8nG f 6^pdv g . (63) 

where the background state is (p-ffir, §<?#)• We thus see the local vanishing of the third term on the 
right hand side when H — > 0. Observe that the kinetic term \q} g deduced from expression 63 (there 
is an extra factor of a half when we read the energy from its Hessian) is consistent with the first and 
second FL equations in the radiation regime, where the densities of gravitational energy and pressure, 

\K\ 2 

are unequivocally identified with pq — pq — J 6 J G . Note however that the first term in 



63 



doesn't 



form part of the effective densities of gravitational energy and pressure in equation 47 and therefore 
doesn't influence the universe's deceleration, instead it is part of the curvature term in the first FL 
equation. 

The calculation of the Hessian is as follows. In terms of conformal variables, a state (g, K) is 
written as 

gab = <fi 4 gY,ab, (64) 

h 

K =tp K Y + -ip g Y . (65) 

,,2 

Where gy is a Yamabe metric of constant scalar curvature Ry = — 6^|- and Ky is a transverse traceless 
tensor with respect to gy. The conformal factor ip must satisfy the Lichnerowicz equation 

Ap+-(p-p 5 )+ l -^p- 7 + 2TTGpp 5 = 0. (66) 

We will take derivatives along a path (g,K)(X) with (g,K)(0) = (-^gn, j;9h), which in turn can be 
seen as a path (gy,Ky,p)(X). Note that tp(0) = 1. Recalling the derivative of the Laplacian ([10 ) 

- (A')/ =< V 2 /, sf>-<Vf,6h>~< V/, dtr g g' >, (67) 

the first derivative at A = of the Lichnerowicz equation is (we are assuming 5^'p = 0) 

Ap' - = 0. (68) 

which shows that p'(0) — identically. Using that fact we get 

V"(0) = ( [ p e dv g y = 6 f p"dv g{0) + [ dv" gy . (69) 



r 6 

's 

Integrating the Lichnerowicz equation and differentiating the integral equation twice gives 
8fc 2 



p"dv g{a) = 2 [ \K' Y \ 2 dv g(a) + 16ttG f 6Wpdv m , (70) 
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from which we get 

6 / y"dv m = / \K'\ 2 dv m + f S^pdv 



9(0)- 



Now let's compute the second term in equation |69| First we note that 

dv„„ 



tr gY 9Y Wy? , ( tr a9' 



+ P£ L )>> 9 



(71) 



(72) 



To compute tr gY g Y we will use the variation formula for the scalar curvature. As the metrics gy are 
Yamabe of scalar curvature —6^- the derivative in A of Ry is zero pointwise, precisely (|10j) 

R' = -A(tr gY g' Y ) + S8g Y - < Ric,g' >= 0. 

Integrating we get 



(73) 



/ < Ric, g' Y > dVg Y 



0. 



(74) 



for all A. Differentiating again at A = we get 



(< Rid,g' Y > + < Ric,g Y > +(R t c ab )(g , Ycd )(g Y c Y(g b Y d ) + (ffic ab )(^ c J(^ c )(^)0d« s( o)(75) 



The Ricci curvature at A = is Ric = 9 gn- Also the functional derivative of Ricci is 
Ric' = l -A L g' - S*(5g') - l -VV{tr g g'). 



(76) 



Observe that from equation 73 we have tr g (mg'(G) = 0. is the Lichnerowicz laplacian and has the 
expression [10] 



A L T ab = V* VT ab + {Ric ac T\ + Ric bc T c a ) - (Rm acbd T cd + Rm bcad T dc ) 



(77) 



Using both facts and also that g' is taken to be transverse we get from equation 75 that 

J^tr gY g Y dv g(0) = 2 Jjg Y \ 2 dv g{0) + ^ J^< 9y, A Lg Y > dv g{0) . (78) 

To compute the Lichnerowicz laplacian we remember that the sectional curvature of 17(0) is — ~, 
therefore 



6k 2 

A L g' Y =V*Vg' Y - —g' Y . 



(79) 



at A = 0. Using the previous equation in equation [78] we get the result of equation [63] after putting 
together equations [72] [70] [69] 
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3.5 The long time ADM limit of the CMC energy: radiation and mass 

gap. 



In this subsection we will introduce assumption (C) and show how, under that assumption, the 
CMC energy converges assymptotically in time to the sum of the ADM masses of the emerging sta- 
tionary solutions plus a radiative term of the radiation in between. The analysis will lead us to argue 



in subsection 3.6 on the validity of the averaging problem in cosmology under assumption (C) and 
asymptotically in time. First we recall the definition of asymptotically flat stationary solution. 



Definition 7 ( Ulf . pg 16) A maximal (k = 0) initial data set (<?, K, N, X) is a stationary asymptot- 
ically flat data iff 

1. (a) goo = -(l-m ) + o(r-% 

(b) 9ij = (1 + 0(r^))S l3 + 0(r" 2 ) ; 

(c) 9^^-e l]k ^x k + 0{r^). 

2. it satisfies the stationary vacuum Einstein equations g = K = 0. 

Now we state the definition of assumption (C). A schematic representation of a space-time (at a 
given time) satisfying assumption (C) is given in figure [I] 

Definition 8 A long time CMC solution satisfies the assumption (C) iff: 

1. (emergence of isolated stationary solutions) after a sufficiently large time there is a finite set 
of pairs of two-spheres (inner and outer) with constant mean curvatures 2/Lq and 2/L(t) re- 
spectively, varying continuously in time (t — 1/H) such that, inside the annulus in between, the 
unsealed flow (g, K, TV, X) decays in the C 1 norm into a stationary solution (go, Kq, Nq, Xq). At 
the outer spheres, \V</> — V0o| < L(t)^t 1 + e • 

2. (the inside of the inner spheres) after a sufficiently long time the volume of the inside of the 
inner spheres grows no faster than t 1 ~ e . 

3. (emergence of the radiative region) after a sufficiently long time the cosmologically normalized 
flow (g,K,N,X), decays uniformly in C , over the exterior region to the outer spheres into the 
flat cone state (gn, — gn , 1/3, 0). 

4- (boundedness of the CMC energy) < — > 0. 

Some remarks are in order. The interior radius Lq is fixed. The exterior radius L(t) grows monoton- 
ically but less than i: lim^oo = 0, in such a way that at cosmological scales the outer spheres 
get smaller and smaller in size. Similarly, the rate at which the solution over the annulus decays into 
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the stationary solution and the rate at which the solution on the exterior region decays into the flat 
cone solution arc left unspecified here. Item four is a global condition that complements the absence 
of explicit decaying rates in assumption (C). In the CMC flow, the interior regions of black holes are 
expected to evolve as tubes of increasingly large size, and therefore increasing volume. Item two gives 
a bound on its grow in the case they form. We want to stress that all these conditions are tenta- 
tive and are not intended to be conjectural. Neither we conjecture a sort of assumption (C) to hold 
generically. The introduction of assumption (C), we believe, provides an starting point in the study 
of the averaging problem directly from the small structure of exact solutions. All these problems are, 
however, difficult problems in the field. Section|4]is an attempt to clarify these issues in pure radiative 
solutions. 

Now let's see how the CMC energy behaves under assumption (C). The second FL equation in 
terms of the CMC energy is 



j ((p + 3p) + ^L)(l + 4>)dv g + E CMC = m 2 J <t>dv g + E CMC , (80) 

where a — — In — k is the logarithmic time. From item four and equation[80]the CMC energy converges 
to the term — 3W 2 J„ <f>dv g with a difference bounded by C/t 1+e . Now let's separate the region of 



integration into the inside of the outer spheres and its outside. Using the Poisson equation 52 we get 



Ecmc = [ < V0,n out > dA+'iU 2 [ 4>dv g + [ ((p+3p)+^)(l+0)cfc 9 +O(t- (1+e >), (81) 
where Hint is the interior of the outer spheres and £l ex t its exterior. Due to item two in assumption (C), 



the second term on the right hand side of equation 81 is an 0(t ( 1+e >). The boundary term approach 



with an error 0(i~( 1+e )) to the sum of the ADM masses of the emerging stationary solutions. We can 
identify the third term on the right hand side of equation [81] is the radiative term because by item 
three <fr — > pointwise on fl ext and the radiation terms from matter and gravitation decouple. Thus 
we get 

Ecmc ~M+1l, (82) 

the total ADM mass plus the radiation energy. This is the same equation as [60] with the additional 
radiative term. A remark has to be said about the radiative term. In an asymptotically flat context 
the ADM energy is a conserved quantity, therefore the radiative contribution to energy measured as 
the difference between the asymptotically Bondi energy and the ADM energy would be a definite 
non-zero amount. In other words there is a definite amount of radiative energy that forms part of the 
ADM energy. In our context, that amount would form part of the radiative term 1Z. Further work is 
needed to show that, indeed there may exist a non vanishing residual radiative energy in the 1Z term. 

We will use the total CMC energy in section [4] to give a rigorous estimation of the gravitational 
energy in the long time for radiative solutions. 
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Figure 1: A schematic representation of the cosmological scale of a universe satisfying assumption (C) 
after sufficiently long time. The emerging stationary solutions are represented with a galactic symbol 
enclosed in a dashed circle representing the outer spheres. The tails coming out from the inside of the 
emerging stationary solutions represent the large tubes developing inside possible black holes. 



3.6 The averaging problem in cosmology. 

We will discuss here the implications of assumption (C) for the averaging problem in cosmology. 
Noting that J\f = we rewrite the second FL equation in the form 

a" _ - f s (4nG(p + 3p) + \K\ 2 )(l + <j>)dv g _ ^ 2 

i{\ + 4>)v 1 + 0' 1 ' 

with 

- ^{ATTG{ P +'ip) + \K\ 2 )(\ + 4>)dv g 



(84) 



m 2 v 

The integrand is the same as in equation [80] therefore we can decompose the integration as we did 



in equation 81 Note that if in equation 80 we write 3Tt 2 J s (j>dv g — — 3tV4>, we get because of item 



four in assumption (C) and the fact that the reduced volume is monotonically decreasing and bounded 
below by V H that (f> = ^gffi + 0(t"( 2+e )) = 0(l/t). This gives the estimation that the factor 1 + 



in the denumerator of equation 83 behaves as 1 + 0(t ). All together gives 



ad^- 3-4nGH- 2 (M ADM +n) { ' ~ + ^ } (85) 

where M.adm is the volume average of the sum of the ADM masses of the emerging stationary 
solutions, and TZ — p ra d + 3prad + PG + Spa where p ra d, Prad and pa, Pg are the volume average of 
the energy and pressure densities of material and gravitational radiation respectively. Equation [85] is 
a differential equation in t, however the estimate on its right hand side is in terms of t — 1/Ti. We 
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thus complement this equation with a differential equation for r as a function of t. From the defining 
equation of r we get the equation 

^ = l + 4> = l-^-t 2 (M AD M + n). (86) 



Equations 85 and 86 are the main equations for the averaging problem under assumption (C) and 
asymptotically in time. We remark that still the Einstein equations have to be used in full, to provide 
an estimation of the radiative term TZ. The next section intends to provide these estimates in the case 
■M-adm — 0, i.e. a purely radiative solution. 

4 Long time smoothing and estimates on the gravitational 

energy: radiation. 

We will use the notation H s for the Sobolev space with s derivatives and iJ| H for the Sobolev 
space where the norms and covariant derivatives are calculated via gn (see [3])- We will prove here 
Theorem [I] The proof is a natural extension of the analysis in [6] . 



Theorem 2 (Expansive smoothing and energy estimates). Let £ be a compact and rigid hyperbolic 
manifold. There is an e > such that the Einstein CMC flow of a cosmologically scaled initial state 
(i.e. with Ti — 1 ) (g,K) with V — V{ n f < e and £ \ < e has the following long time properties (take 
t = 

1. The limit Umt—>oofiQa is finite and greater than zero. 

2. there are ni > such that lim^oo jj^jyr-Qi < oo for i > 1. 

3. for given 7 > 0, / t °° S ^ dVa du > C^( 2 +^. 

4- \K\ 2 ^ Ci -4 pointwise (not volume averaged). 

In particular the cosmologically scaled flow of a H l x H 1 ^ 1 state (for any i > \) as in the hypothesis 
above converges in H l x H 1 ^ 1 to the canonical flat cone state (g,K) = (gui^da)- 

Proof of theorem [l] We start by recalling a result from [B] that will be useful to prove items 1 and 
2 in theorem [TJ 

Lemma 1 [T] Let £ be a compact and rigid hyperbolic manifold. There are C and £q such that if a 
cosmologically normalized CMC state (g,K), where g is harmonic with respect to gn, is e-close to 
{gn , —gu) in the Hg H X Hg H topology, with e < £o then there is a constant C (dependent on eo) such 
that 

< (\\g-g H f H3 +\\K + g H \\ 2 H2 )<C£ 1 . (87) 

SH 9H 
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We get therefore the elliptic estimate for the Newtonian potential <f> = N = — 1 from the lapse 
equation 

\\N\\h, b < C\\K\\ H 2 h \\K\\ L 2 h < CEi, (88) 



and 



\N\\m < C\\K\\ H 2 \\K\\ H i < C£ x . (89) 



To extract conclusions on the decay of the Sobolev norms of the cosmologically normalized states 
we will make use of the fact proved in [3] that under the conditions of the last lemma, eo & n d £i-\ 
controls the difference of the states in H gH x H^" 1 with respect to the background state (gH,—gii) 
states at zero, i.e. the derivatives tend to zero in L 2 gH as eo and £i tend to zero. 

Item 1. The Gauss equation gives the following inequality for the evolution of the first order cosmo- 
logically normalized Bel-Robinson energy ([B]) 

^<-2£ 1 + C£?. (90) 

OCT 

with c a constant greater than zero. It follows therefore that £\ decays faster than the solution x(a) 
to the following ordinary differential equation and same initial condition 

x' = -2x + cx%. (91) 

This is a Bernoulli type of equation that can be solved by making the change of variables v — 
which gives the differential equation 

«' = «-|, (92) 
having the solution v — ^ + Ae® . This implies that 

x = *(q )e-'(—°> (93) 

(|( e -(a- ffo )_l)a;(o-o)i + l) 2 ' 



which results in the following decay of £ 



- (c (e -(a-ao)_ 1)^(^)1 +1)2' 

Observe that if Co is big enough then we get the bound 

* < ^ 4 2i -°\ ( 95 ) 

Now we prove item 1 in theorem [T] From the Gauss equation and Lemma [T] and the above estimate 
for £ i we get an evolution equation for Qq of the form 

^ = -200 + ^), (96) 
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where h(a) is a function which is bounded in absolute value by 

|%)|<4(a )e- 3 ( ff -'»), (97) 
Therefore we get the following expression for Qq 

Qo = e- 2(CT - CTo) (QoM + e- 2ao f h(u)e 2u du), (98) 

Clearly the integral in h has a limit when a — > oo. If the term in parenthesis on the right hand side 
has limit different than zero then we are done as then 

lim > 0. (99) 

Let's see that the limit cannot be zero. If that happens then we have for all a 

/•OO 

Qo(o-) = -e- 2a \ h(u)e 2u du. (100) 



The integral is negative for all a (Qq is positive) and goes to zero as a — ► oo. Then there is a diverging 
sequence {(Ji\ such that for all a > Oi we have 

/•oo 

h(u)e 2u du <- h(u)e 2u du (101) 

J (Ti 

making then 

Qo(<r)<0o(o-»)e~ 2(<r ~ ffl) , (102) 

for all cr > ai. Using again the Gauss equation, Lemma [l] and the estimate above we get an evolution 
equation for Qo(a) of the same form as in equation [96] with h instead bounded in absolute value by 
C£i(ai)^ Qo(<Ji)e~ 3< - r7 ~ (7i \ It thus gives an expression for Qo of the form 

Q (a) = Q Q {ai)e- 2 ^\l + f k ^-^du). (103) 
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To see that linio-^oo Qoe 2a > we note the following bound for the integral term in the equation 
above 

v _ 2ai r^^^^^ (104) 

which tends to zero as ct; — > oo. This is a contradiction, thus the limit must be positive. 

Item 2. Now we prove item 2. By induction we will be able to get an equation for £i(a) of the form 

Q^ = -(2 + /i'(a))Q i + /i((7)Q?, (105) 
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where h'(a) and h(a) are functions bounded in absolute value by C'er™ e a and Ca n e a for some 

„ i 

C",C and n',n constants. It follows after making the change of variable v — Q? that Qi can be 
bounded by an expression of the form 

Qi < C(T 2 (" +1 )e- 2CT , (106) 

for some constant C . 

Lemma 2 Suppose that a solution to the CMC flow (<?, K) has 

Qj(o-) < Cj(j nj e^ 2a , (107) 



for j — 0,...,i > 1, then Qi+i satisfies an equation of the form 105 and therefore satisfies an 
asymptotic of the form \T07\ for j = i + 1. 

Proof: We start with the differential inequality for Q l . Make /? = ^p. Then Qi(k) = A^^Q^fc), 
and therefore 

A useful trick for the calculations that follow is to write 

where the f3 inside the derivative on the right hand side is taken constant equal to its value at the 
time of differentiation. Thus we are calculating the fc-derivative of the cosmologically scaled solution 
at k = —3. Putting all together we get 

— = ( 2i + i)g 4 + 3 d{(jk) . (110) 

We are going to study the derivatives of perturbation of the canonical flat cone state (<?# , — gn) 
at k — —3. From the Gauss equation we have 

dQ(i 
~dk 

therefore 



-3 / NQ {l]abTT n ab dv g - / 2N(E$J (r)aTb + B§J* (i)aTb )dv g , (HI) 



E 



3/3 2l+1 ^ = -9 ^ JVQ WoKff fl ab dv g - £ 6A>(£?" J ( . )af 6 + BtJ* {€)aib )dv~ g . (112) 

We will say that a term is an O(a) if it can be bounded in absolute value by a term of the form 
Ccr"e~ <T for some natural number n. Let's start by analyzing the first term on the right hand side of 



equation |111| Making 

n ah = n ab + ^(g a6 + T a T b ), N = N-^, (113) 
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we get 

-9 / NQ abTT fl b dvg - 3Q t - 9 f NQdvg. (114) 
Using Lemma [l] and the estimate on £y above we get the term 

-3Q 4 + 0(ct)Q 4 . (115) 



Now we estimate the second term in equation 111 and therefore we need estimates of J and J*. We 



will do the calculations only for J, those for J* proceed in exactly the same way. We note first the 
following inductive formula for J 

J(Wi) abc = n de V e W (l _ 1)dafcc - IW (l)dabc T d + T * Rm * W^i+Vr J(W*) abc (116) 

where the * is some tensorial multiplication whose particular form is not important to our purposes. 
We can write the formula above symbolically as 

J(W<) = n*VWi_! - fW t - * T - |J(W,_0 + T * Rm * Wj-i+VrJCWj-i). (117) 

Now, inducting the fifth term on the first, second, third and fourth gives the following terms respec- 
tively 

i=i— l 

i. ^(n*vWi-i-i) (us) 

j—'l— 1 

2 - E V t(^* T * W ^'), (H9) 

»— 2 

3- E V t(tT J ( W *-«+i)))> (12°) 
i-l 

4. ^V 3 r (T*Rm*W j -i_ i ). (121) 

~ i 

The only terms that are not going to count as O(a) or 0(a)Q? are those coming from the expression 
2 and when the Vt derivative applies only to the Wi_j giving 

yiW,*T (122) 

When we take into accound this and a similar term aroused from a formula for J* and plug them into 
equation |111| we get a contribution of the form 

-2iQ, (123) 
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As said above and as we will explain in a moment all other terms are going to count as O(a) or 



0(cr)Q? therefore we would get, putting equations 110 115 and the last estimate together we get 

d ^ r = -{2 + 0{a))Q i + 0{a)Qh (124) 

as we wanted in the induction. To discuss the other terms then we start by recalling some propositions 
from [5] restated in a different form for convenience of the article. 

Lemma 3 Let (g,K) be a CMC flow on a rigid hyperbolic manifold S. Suppose that the initial cos- 
mological state is e-close to the standard flat cone state (<?#, ~9h) as i- n Lemma^then (all derivatives 
below are taken at k — —3 ) 

1. ||V^n|| ffJ , i> 1, j = 0,1,2, (125) 

SH 

are controlled by Si+j-±. 

2. (HVWOHj* + HWilUg ) < C(||W i+1 |U + ||W,||m + || J(Wi)|Uj ) (126) 
i > 0. 

Lemma 4 V^J(W;) has an expression of the form 

v^j(Wi) = xxv^n)™ 1 * • • • * (v™ s n)™* * n'*vWi + (127) 
+ ^(v^n)" 1 * • • • * (v^ s n)"'- * n r *v^(T * Rm * w £ ) (128) 

where the first sum is among the set k < i+h — 1, mi > . . . > m s > 1 and ■ 7ij(l+TOj)+Z+fe = 
while the second is among the set rh\ > .. . > m s > 1 and n.j(l + rhj) + k + l + q = i + h — 1. 

Now we prove the following Lemma. 

Lemma 5 Let (g, K) be a CMC solution. Suppose for a given value of i there are rij and Ci such 
that Si < C l o- n *e- 2a = O(cr) then 



1. there are n\ and C[ such that ||J(W i )||? 2 < C' l a n 'e^ 2a = O(a). 

SH 

2. There are n' (J - and C' VJ such that ||(V 3 T J(Wi_ j ))~|||a < C' VJ a n '^ e - 2a = O(a) for j < i. 

Proof: 1. Proceed by induction in i. Observe that all the factors involving II and its time derivatives 
(with h = 0) are controlled by Si in H 2 by Lemma 3 The norms ||(VWfc)~|| L 2 are 



in formula 



127 



controlled using inequality |138| The second kind of terms in equation |128| are controlled as follows. 
The factors involving II and its time derivatives are controlled again in H 2 H by Si. The other factors 
can be seen as 

(Vf,(T * Rm * W s ))~ = E qi+q2+q3=q (^ q TTr * (V??Rm)~ * (V^W S )~, (129) 
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= ]T (V£T)~ * (W 92 ) * (W q3+ ~ k ), (130) 

91+92+93=9 

with q < i — 1. Now Sobolev embeddings give 

l|W, 2 *W £+ , \\„ <C(\\W q3 \\ Hl J\W~ k+Q \\ H i J, (131) 



2 

9H 



where the factors on the right axe controlled by Lemma [5] The factors CVtT)~ are controlled in H : 
by Lemma [3] Finally the proof of part 2. is the same as above after using formulas |127[ [T28] □ 



The terms in 2, 3, 4 on the induction formula for J other than the ones already considered in 
I - 1 

equation 123 are easily seen to be bounded by O(a) or 0(a)Q^ by the same kind of arguments as in 

Lemma [5] To bound the terms in 1. in the same way we need the following form of V^VWj. 

v 3 T vw, = ^(v™ 1 ^" 1 (v™ s n)"*- * n'*vw t + (132) 

+ ^(V™ 1 !!)™ 1 * • • • * (V™ S 1T)"- * n [ *V£(T * Rm * W s ), (133) 

where the first sum is among the set mi > . . . > m s > 1 and • n j(l + m j) + 1 + k — i + j, while the 
second is among the set fh\ > . . . > m s > 1 and ■ rij(l + m^) + fc+Z + gr = i + j — 1, which can be 
easily proved by induction by using equation 

V T VWj =VW l+1 + n*VWj + T * Rm * W 4 . (134) 

This finishes the induction in Lemma |2] □. 
Items 3 and 4. The estimate from above in item 4 comes from Lemma [TJ The item 3 or the estimate 
from below is more involved, the argument is as follows. 

Lemma 6 For any e > there is a ball Br gH _ gH \(8) of cosmologically scaled states in H 3 x H 2 such 
that 

\\N-\\\ L -<e, (135) 



and 



1 



A-kGHEcmc > ^ J s \K\ 2 dv g . (136) 

We can prove item 3 by making use of the Lemma [6] First, the derivative of the reduced volume 
V = H^V in logarithmic time is 

^ = -3 J^N\K\ 2 dv g . (137) 

If we integrate it from a to oo and use lemma [6] above we get the following inequality 

— J— / \K\ 2 dv g < 4-kGHEcmc = 3 / ( / N\k\ 2 dv g )da < (1 + e) / ( / |l| 2 d« g )d<T. (138) 
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Making U — f°°{f s \K\ 2 dvg)da the inequality 



138 



is written as 



U' > -(4 + e)(l + e)C/, (139) 
which after integration gives the left hand side inequality in item 3. 

Proof of Lemma [(i] First we note that the estimate for N — | is deduced from Lemma [TJ For 
the second estimate it may be deduced from the calculation of the Hessian of the energy that we did 
before, however we will follow a direct estimate from the Lichncrowicz equation. We argue as follows. 
Say g — </> 4 <?y where gy is the unique metric in the conformal class of g having scalar curvature —6. 
Then cj> satisfies 

_A0+^ 5 -0) = V 3 |#|y- (140) 
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The maximum principle gives <f) > 1. Making = 0—1 rewrite equation 140 



as 



-A0+^(0 3 + 2 + 0+l)0=^f. (141) 
At the point where 4> or (j) is maximum we have 

" 12 03 + 02 + 0+1 <^2". ( 142 ) 
which gives if ||-ftr||i=o is small 

ll^llioo 

U\\l°° < • (143) 

Wl1 - 12 - \\K\\l T 1 ' 

Also note that 

-o-(E) < -6V^f, (144) 
which gives 

< / (0 6 - l)dv Y <V~V H . (145) 

s 



Writing 6 - 1 = (0 - l)(</» 5 + 4 + 3 + 2 + + 1) we get 

6 / {4>-l)dv gY <V-V H . (146) 
Integrating equation |140| we get 



/ (0 5 - <t>)dv gY = / 0- d |tf|^ Sy . (147) 
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Under the assumptions we have and using equation |143| we can get from equation |147| above the 
inequality 

6(4 + e) f (4>-l)dv gY > f <j>- 2 \k\ 2 Y dv gY = / \K\ 2 dv g . (148) 



which together with equation 81 gives the inequality 



(4 + e)(V-V H )> / \K\ 2 dv g . (149) 
Js 

as desired. □. 

This finishes theorem [l] □. 

5 States of arbitrarily large gravitational energy. 

We will construct a one parameter family of states (g\,K\) such that 

1. k\ — ko fixed, 

2. Vol gx -^a^oo oo and \\K\\\ L 2^ -^a^co oo, 

3. The "big-bang" family of states, i.e. the volume-one normalized family of states above has 

-fc A ^oo, (150) 

Vo/ flx (E) = l, (151) 

lim \\K X \\ L 2 =oo. (152) 

As has been argued above, these states represents a one parameter family of states with arbitrarily 
large gravitational energy. The construction is as follows. Pick the hyperbolic metric g^ and a non 
zero transverse traceless tensor K with respect to it. According to the conformal method it is possible 
to find a solution to the constraint of the form {g\, K\) = {ip gu, \ 2 ip~ 2 K—(pgH) (the mean curvature 
being k = ko = —3 and one parameter family of states as above with arbitrary ko can be obtained by 
scaling), by solving the elliptic equation 

A^=-^-y|^lL^ 7 +^ 5 - (153) 
Now we prove items 2. Multiplying equation |153| by <p and integrating we get 

y ^ \K\* BB <p-*dv ga = £ |V^| 2 + - <f)dv m . (154) 
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Note that the left hand side is | |jX;>J|^2 . If the left hand side doesn't diverge as A — > oo then the 
right hand side remains bounded in particular the H\ H norm of p remains bounded. Pick an open set 
f2 where |-K"| Sff > e > 0. Then as ip is bounded in H 1 we have Vol{x S Q,/p(x) < n} — > Vol(Q) as 
n — > oo uniformly in A. Then for some n we have Vol{x G fl/ip(x) < n} > ^"^^ uniformly in A, and 
so the left hand side is bigger than e 2 Vol(Q) which diverges when A — > oo which is a contradiction. 
This proves item 2, to prove item 3 we argue as follows. The L 2 norm of K\ of the volume one states 
are 

A 4 / s |iT|V^ 3H = J s |Vyf + |(y 6 -y 2 )<fo flH . 

We have that an upper bound on the left hand side in the last equation implies an upper bound for 
the H 1 norm of ip, for if not we have J s ip e dv gH — > oo which would make the numerator of the right 
hand side diverging in A, but we know J s ip 6 dv gH diverges which is a contradiction. 



6 Summary and open questions. 

We have introduced the notion of general JC = — 1 cosmological model as a formal definition al- 
lowing to study cosmological notions in arbitrary solutions of the Einstein equations. This gave us 
a framework to study general cosmological solutions in a cosmological language. The approach may 
be applicable to models other than general JC = — 1 cosmological models, i.e. models with different 
spatial topologies. Thinking on the averaging problem in cosmology we have defined volume-averaged 
cosmological parameters and an averaging map: a correspondence between arbitrary solutions and ho- 
mogeneous and isotropic Lorentzian spaces. Those concepts allowed us to give a precise mathematical 
formulation of the averaging problem in cosmology. In another section and aiming at the start of a 
rigorous analysis of cosmological evolution from the solutions at the natural scale, (i.e. including the 
small scale), we have introduced assumption (C) which precisely describe a certain class of solutions. 
Those solutions are divided into two main subclasses: radiative and mass gap. We have given a de- 
tailed description of the full structure of the radiative solutions. We have also analyzed the averaging 
problem in cosmology in precise quantitative terms for mass gap solutions. The attemp may be con- 
sidered as a first step towards the ideal goal of attacking the averaging problem in cosmology directly 
from the solutions at the small scale. Finally we constructed initial "big-bang" states of arbitrarily 
large gravitational energy, showing that, apriori there is no mathematical restriction to assume the 
gravitational energy to be low at the beginnings of time. 

There are several questions and avenues of research left open in the present article, of varying 
difficulty however. For instance one may want to see in action the formalism of general cosmological 
models in cosmological solutions with Cauchy surfaces of non hyperbolic topology. Also and perhaps 
more important is to obtain rigorous results that may support or not assumption (C). Any rigorous 
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result of the sort would put the study of the averaging problem in cosmology from the small scale on 
a firm basis. Analyzing the validity of assumption C from the Einstein equations is a very difficult 
problem. A central point is to study the spatial asymptotic of stationary solutions that may emerge in 
time. Is an emerging stationary solution necessarily spatially asymptotically flat in the long time?. If 
the answer is affirmative one may be in a better position to prove the a priori estimates in assumption 
(C). The answer may instead be negative and that would open a new avenue of research. Finally the 
analysis of the validity of the averaging problem in cosmology from assumption C was only asymptotic 
in time, and therefore of non obvious applicability. An interesting question is to study the validity of 
the analysis but in finite times. 
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